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Abstract 

A nonlinear equation describing curved stationary flames with arbitrary gas 
expansion 9 = pfuei/Pburnt) subject to the Landau-Darrieus instability, is obtained 
in a closed form without an assumption of weak nonlinearity. It is proved that 
in the scope of the asymptotic expansion for 9^1, the new equation gives the 
true solution to the problem of stationary flame propagation with the accuracy 
of the sixth order in — 1. In particular, it reproduces the stationary version of 
the well-known Sivashinsky equation at the second order corresponding to the 
approximation of zero vorticity production. At higher orders, the new equation 
describes influence of the vorticity drift behind the flame front on the front 
structure. Its asymptotic expansion is carried out explicitly, and the resulting 
equation is solved analytically at the third order. For arbitrary values of 9, the 
highly nonlinear regime of fast flow burning is investigated, for which case a large 
flame velocity expansion of the nonlinear equation is proposed. 



1 Introduction 

Curved flame propagation is one of the most important and difficult issues in the 
combustion theory. Despite considerable efforts, its closed theoretical description is still 
lacking. Perhaps the main reason underlying complexity of the problem is the intrinsic 
instability of zero-thickness flames, the Landau-Darrieus (LD) instability |l|, @]. In view 
of this, evolution of the flame front cannot be prescribed in advance. Instead, it should 
be determined in the course of joint analysis of the flow dynamics outside the flame 
front and the heat conduction - species diffusion processes inside. In general, nonlinear 
interaction of different perturbation modes under the smoothing influence of thermal 
conduction leads to the formation of a steady curved flame front configuration with the 
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curvature radius of the order lOLj, where Lf is the flame front thickness. This estimate 
for the curvature radius can be obtained from the hnear theory of the LD-instabihty 
0, where it corresponds to the cutoff wavelength for the front perturbations; it is also 
confirmed by 2D numerical experiments on the flame dynamics 

In essential, difficulties encountered in investigation of the nonlinear development 
of the LD-instability are twofold: 

1) Perturbation analysis of the nonlinearity of flame dynamics with arbitrary gas 
expansion is generally inadequate. In particular, it is completely irrelevant to the 
problem of formation of the stationary flame configurations. 

2) Finite vorticity production in the flame implies that the flow dynamics down- 
stream, and therefore the flame front dynamics itself, is essentially nonlocal. The latter 
means that the non-locality of equations governing flame propagation is more complex 
than that encountered in the linear theory and described by the Landau-Darrieus op- 
erator. 

Only in the case of small gas expansion can the problem be treated both perturba- 
tively and locally, since then the amplitudes of perturbations remain small compared 
to their wavelengths at all stages of development of the LD-instability [|], 0, |H, and 
the flow is potential both up- and downstream in the lowest order in ^ — 1, where 6 is 
the ratio of the fuel density and the density of the burnt matter. In this approxima- 
tion, the nonlinear evolution of the front perturbations is described by the well-known 
Sivashinsky equation 

The nonlinear dynamics of flames with finite gas expansion has been the subject 
of a number of papers some of which are briefly considered here as the characteristic 
examples of dealing with the difficulties mentioned above. 

To render the problem tractable, certain simplifying assumptions has been intro- 
duced by various authors in attempt to weaken or even get rid of one of the two above 
features inherent to the nonlinear flame dynamics. In Ref. 0, the vorticity production 
in the flame is completely neglected [see the point 2) above]. The mass conservation 
and the constant normal flame velocity are taken as the conditions governing flame 
dynamics. The problem is reduced thereby to the well-known electro dynamic problem 
of determining the single layer potential with constant charge distribution proportional 
to the gas expansion. However, neglecting the vorticity downstream breaks the conti- 
nuity of tangential velocity components as well as the constant jump of pressure across 
the flame. This model thus violates the basic conservation laws to be satisfled across 
the flame front. Note that the number of neglected degrees of freedom (three in the 
3-dimensional case) just corresponds to the number of broken conservation laws (two 
tangential velocity components and the scalar pressure). 

In contrast, in Ref. an attempt is made to take into account the vorticity 



production in stationary flames, under the assumption of weak nonlinearity [see the 
point 1) above]. From the mathematical point of view, however, the assumptions of 
weak nonlinearity and stationarity contradict each other. Using them simultaneously 
turns out to be inconsistent except for the case of small gas expansion. Indeed, let 
us consider a weakly curved flame front propagating in z-direction with unit normal 
velocity^ with respect to an initially uniform fuel; the transverse coordinates will be 
denoted x. It is not difficult to show that in this case the flow equations up- and 
downstream together with the conservation laws at the flame front imply the following 
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relation between the flame front position, z = /(x, t) — t, and ^-component of the fuel 
velocity v = v(x, z, t) just ahead of the front, m„ = v^(x, /(x, t) — 0,t), 

+ ^ = ^*^' 

where $ denotes the Landau-Darrieus operator defined by 

/ + 00 
rfk|k|/kexp(2kx), 
-oo 

/k being the Fourier transform of /. (It is assumed here that \{z = — oo) = in the 
laboratory frame of reference. Details can be found, e.g., in Ref. |TT|) On the other 
hand, to the leading order in the front slope df /dx, the condition of unit normal flame 
velocity with respect to the fuel gives 

In the regime of steady flame propagation, 

dt ' dt 

where W ~ (9//9x)^ is the flame velocity increase due to the front curvature. Substi- 
tuting this into relations (P, (||) gives 

Since the right hand side of this equation is linear^in the slope, while the left hand side 
is only quadratic, nonlinearity can be considered small only if ^ — 1 is small, in which 
case one has df /dyi ~ ^ — 1, ~ (^ — 1)^. Thus, for arbitrary 9, the weak nonlinearity 
approach to the stationary flames, advocated in [|10] , turns out to be self-contradictory. 



From the physical point of view, this means that a weakly curved flame with finite gas 
expansion cannot be stationary. Instead, equations (|l|), (|) should be considered in this 
case as determining evolution of the small perturbations in time. Indeed, substituting 
Eq. (H) into Eq. (|l]), we obtain, in the linear approximation, 

which is nothing but the well-known LD dispersion relation for the perturbation growth 
rate, written in the coordinate space 0]. 

In practice, the discrepancy in the orders of magnitude, expressed by Eq. @, 
shows itself as the impossibility to correctly develop perturbation expansion in powers 
of the slope. For instance, the jump of the pressure field across the flame front, found 
in Ref. [|10], Eq. (33), turns out to be velocity-dependent already for zero thickness 
flames, while according to the well-known exact result [0 this jump is constant along 
the flame front, and is equal to 1 — 6. 

^See the Appendix. 
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It should be stressed that the above-mentioned inconsistency resides in the equa- 
tions of motion themselves; it is inherent, therefore, to the flame propagation in tubes 
as well as to unbounded flames. Indeed, even in the case of very narrow tubes, numer- 
ical experiments on flames with 6 = 6^8 give values of about 1.5 2.0 for the slope 



Finally, we mention Ref. [|13[, where a non- stationary equation describing the inter- 
action of perturbations at the early stage of development of LD instability is obtained 
at the second order of nonlinearity, under a certain model assumption concerning the 
structure of downstream flow. Namely, it is assumed that there exists a local rela- 
tion between the pressure field and the potential mode of the velocity downstream 
(expressed by the Bernoulli equation). However, as we mentioned in 2), the fiow dy- 
namics downstream is essentially nonlocal, in particular, the pressure field is expressed 
through the velocity field by an integral relation, its kernel being an appropriate Green 
function of the Laplace operator. From the work |]13|, one can see once again that 



the perturbative treatment of nonlinearities is not applicable to the stationary fiames 
with arbitrary 6 : since, at the first order, the time derivatives of the front position 
are linearly related to its space derivatives through Eq. (^), stationary form of the 
non-stationary Zhdanov- Trubnikov equation depends on the way this relation is used 
before the time derivatives are omitted. Transition to the stationary regime in this 
equation is therefore ambiguous. 

In view of what has been said, we arrive at the conclusion that the stationary fiames 
can only be treated consistently in the framework of the small 6 — 1 expansion. Thus, 
the problem of describing the stationary fiames, and more generally, non- stationary 
fiames in the fully developed nonlinear regime, can be formulated as the problem of 
deriving an expansion of unknown exact equation, governing fiame djTiamics, in powers 
of e - 1. 

The purpose of the present paper is to show that in the stationary case, such an 
expansion can indeed be developed beyond the second order, the order of validity of 
the Sivashinsky equation. We found it convenient to employ the model assumption of 
Ref. in our analysis, which allows one to obtain an equation describing fiames with 
arbitrary gas expansion in a closed form, without the assumption of weak nonlinearity. 
This will be shown in Sec. |3.1| on the basis of simple geometric considerations. That the 
equation obtained actually provides the true expansion of (unknown) exact equation 
with the accuracy of the sixth(!) order in 6* — 1 is proved in Sec. |3.2| . This expansion 
is carried out explicitly in Sec. [4.1| . The third order equation turns out to have the 
functional structure analogous to the Sivashinsky equation, and therefore can be solved 
analytically. This solution is found in Sec. [4.2| . Finally, the model equation obtained 
in Sec. |3]1| is invoked in Sec. ^ for investigation of the highly nonlinear regime of fast 
fiow burning, where a large fiame velocity expansion of this equation is proposed. The 
results of the present work are summarized in Sec. ^. Some mathematical results used 
in the text are derived in the Appendix. 



2 Flow equations and conservation laws 

Let us consider a 2D fiame propagating in an initially uniform premixed fiuid. Since 
our main concern is the infiuence of finite gas expansion on the fiame structure, we 
will assume in what follows that the fiuid is an ideal gas with constant thermal con- 
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duction and constant specific heat. The ratio of the heat and mass diffusivities (the 
Lewis number) is taken to be unity. For definiteness, we consider flame propagation 
in a tube of arbitrary width with ideally shp adiabatic walls. Under the assumption 
that development of the LD instability ends up with the formation of a steady flame 
configuration, we go over to the reference frame of the stationary flame. Let the Carte- 
sian coordinates (x, z) be chosen so that 2;-axis is parallel to the tube walls, z = — oo 
being in the fresh fuel. It will be convenient to introduce the following dimensionless 
variables 

where Uj is the velocity of a plane flame front, P/ is the initial pressure in the fuel far 
ahead of the flame front, and R is some characteristic length of the problem (e.g., the 
cut-off wavelength). The fluid density will be normalized on the fuel density p_. As 
always, we assume that the process of flame propagation is nearly isobaric. Then the 
velocity and pressure flelds obey the following equations in the bulk 

du dw ^ , , 

du du 1 dU . , 

dw ^ dw 1 dU 
d^ drj p drj 

The above flow equations are complemented by the following conservation laws to 
be satisfled across the flame front 

u+ - 9u- - f'{w+ - OwJ) = 0, (8) 
W+-W- + f{u+ - U-) = 0, (9) 

n+-n_ = -(^^-i), (10) 

where the flame front position is given by ^ = f{i]), the subscripts and "— " mean 
that the corresponding quantity should be evaluated at ,^ = /(//) -|- and ^ = /(??)— 0, 
respectively, and 

df 



Finally, the normal velocity of the fuel at the flame front is unity, which is expressed 
in the form of the evolution equation 

u_ - fw_ = N. (11) 

Equations (|8|)- (p!ID are written for zero thickness flames. These are of primary interest. 



since in the majority of cases the thickness of the flame front is small compared with 
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the fluid-dynamical scale of the problem. Formal generalization to the case of small 
but nonzero thickness is straightforward and will be done in Sec. |3.3| . 

In principle, the system of equations (|^) - (|Tl|) completely determines stationary 
flame configuration as well as the flows upstream and downstream. To derive an 
equation for the flame front function f{T]), one has to find solution of the bulk equations 
(H) ^ ® for the fuel and combustion products with appropriate boundary conditions 
on the walls, subject to the conservation laws - (p!OD which are a kind of boundary 
conditions on the flame front. The requirement that the obtained solution satisfies 
Eq. (pUD then gives an equation for the function /(r/) itself. 

As to the flow upstream, the corresponding solution is readily obtained: since the 
flow is potential at = — oo {u = V, w = 0, where V is the velocity of the flame in 
the rest frame of reference of the fuel), it is potential for every C, < f{fi) in view of the 
Helmholtz theorem thus 



u 



2 

where the Hilbert operator 



-t-oc 

V + u = V+ j dk Ukexp{\k\^ + ikr]), (12) 



w = Hit, (13) 



n + ^(^^^ + w^^) = consti, (14) 



1 T.. /(c) 



{Hf){r^) = -p.v. / d( 

TV J 



-oo 



"p.v." denoting the principal value. Although the relation w = Hu between the velocity 
components upstream is nonlocal, it is expressed in terms of the transverse coordinate 
7] only. 

Things become more complicated downstream. There, no relation exists for the 
variables m, w, 11, which can be expressed in terms of rj alone, since the ^-dependence of 
these variables is unknown because of the presence of vorticity produced by the curved 



flame. Nevertheless, we will assume following Ref. [13| that a potential mode v. 



[w. 



pi 



Up) can be extracted from the downstream velocity v = (w,m), such that the 



following Bernoulli-type relation holds between Vp and 11 

1 



n + ^(wJ + «^') =const2. (15) 



Then Eqs. (^) - (^ can be rewritten as 



du dw ^ dup dwp ^ 
dC, drj ' dC, drj 



du du dup dup 
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dw dw dwp dwp 
The general solution for the potential mode can be written analogously to Eqs. (|T^, 



(0) 



Up = 9V + Up = 9V+ J dk Ukexp{-\k\^ + ikT]), (19) 



Wp = -Hup. (20) 

The model relation (pTS]) does not uniquely define the potential mode, since a constant 
term in the ^-component of the velocity can be assigned either to the potential mode 
or to the vorticity mode. As will be shown in Sec. |3.2| , the vorticity produced in a 
fiame is only of the fourth order in (^^ — 1) as — 1. The choice 9V of the constant 
term in Eq. (^^ is fixed therefore up to the second^ order by the requirement that the 
mass fiow at ^ = +oo, where u = Up = const, = 0, equals that at ^ = —oo. 

Now, using the continuity equations (|TBp for v and Vp, Eq. ([T^) can be written in 
the form in which ^-dependence is implicit 

dw du dwp dup 

w- Up— \-Wp^— = 0. (21) 

or] or] or] or] 

It will be shown in the next section, the above equations for the potential mode up- 
stream and downstream, equation (pTD, and the conservation laws at the fiame front 
constitute the system of equations sufficient to derive an equation for the function /(r/) 
in a closed form. 



3 Nonlinear equation for the flame front 
3.1 Derivation 

It will be shown presently that the set of equations (H) - ([TT|), (|13D - (llSD, (^0|), and 



(pl|) can be transformed into one equation for the function f{r]). The fact that this 
set is written in a form that does not explicitly operates with the ^-dependence of the 
fiow variables makes it unnecessary to follow the program outlined in the preceding 
section to obtain an equation for fij]). Namely, the specific structure of the up- and 
downstream fiows in the bulk is now irrelevant. In particular, knowledge of the ^- 
dependence of the velocity and pressure fields turns out to be superfiuous. Roughly 
speaking, the ^-dependence of a function F(^, r]) describing the shape of the fiame front 
is known in advance, since the equation F{^,r]) = can always be brought into the 
form C,—f{r]) = (with / many- valued, in general). Determination of the //-dependence 
alone of the functions involved is therefore sufficient for the purpose of description of 
the fiame front structure. 



■^One might think that since the vorticity is of the fourth order in 6* — 1, the flow is to be potential 
up to the third order. One should remember, however, that the expression of the velocity field through 



the vorticity field is nonlocal. Specifically, we will see in Sec. 3.2 that the vorticity produced in the 
flame can be expressed as the ry-derivative of a certain function of the velocity, while differentiation 
along the fiame front brings in an extra power of 6 — 1. 
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In what follows, it will be convenient to introduce separate designations for the up- 
and downstream velocity and pressure fields. Namely, they will be distinguished by the 
supersciipts " +" and " — ", respectively. Then, setting = firf), equations ([l3|) , (p!^ , 
(p^), (pOD, and (|21|), together with the conservation laws (|D - (|10|) and the evolution 
equation ([TT|), can be rewritten identically as follows 



— 9u - 



f'iw 



n+ 



25 {(""p 



U 



U 



— w 



- du+ 



U 



. dwp 
P dr) 



+ -ew-) = 





— U^) = 





n+-n = 




--f'w- = 




nj- = 


Hu- 


+ iw-r} = 


consti 


+ K)n = 


const2 




-m; 


1 -1- duf 

+ = 

P or] 






«=/(r?) 



Suppose we have found a solution / = firf), v~ = v~(^,77), v"*" = ^^{i,rf), etc. of 
the set of equations in the large brackets in (*). Then, in particular, these equations 
are satisfied for ^ = f{ri). On the other hand, since no operation involving ^ appears in 
these equations,^ the function f{ri) is one and the same for all solutions. Furthermore, 
since f{ri) is ^-independent, it is convenient to work with the particular solution in 
which all the other functions are also ^-independent, and to omit the large brackets in 
(*). Therefore, we replace the above set of equations by the following 



V 



dr] 





9v- - f{uj+ - Ouj-) 


= 


(22) 




-uj- + f\v+ -V-) 


= 


(23) 




TI+ — 




(24) 




V- - fuj- 


= N 


(25) 






= Hv- 


(26) 


7!-- 




= const 1 


(27) 


71+ - 




= const 2 


(28) 




UJ+ 


= -Hd; 


(29) 


j^dv+ 
dr] 


, duo+ , dv+ 

V+ + UJ+ ^ 

^ dr] ^ dr] 


= 0, 


(30) 



where v,u, and vr are the ^-independent counterparts of the flow variables u,w, and 
n, respectively, and 



V 



V 



V 



V,. 



ev + v: 



(31) 



The fact that now the function f{r]) does not enter the arguments of these variables 
allows us to avoid expanding them in powers of / (employed, e.g., in Refs. |10])- 
In fact, such an expansion is irrelevant to the issue whatever regime (stationary or 

^No such operation can appear in the boundary conditions to these equations neither. Otherwise, 
steady flame propagation would be impossible. 
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not) is considered, since all the equations governing flame propagation are invariant 
with respect to the space translations, and therefore, all terms containing powers of 
undifferentiated / should appear in invariant combinations in the final equation for /. 
On the other hand, in view of the above-mentioned translation invariance, the function 
/ itself does not need to be small even if the front is only weakly curved (e.g., at the 
early stage of development of the LD- instability). We thus see that the /-dependence 
of the flow variables through their arguments must eventually cancel in some way. 

Let us now turn to the derivation of the equation for the function f{r]). 

From the geometric point of view, Eqs. (0), (|2^), p8| ) determine the amplitude fl 
of the complex function 

v+ + = ne^'f', (32) 

while Eq. ( ^OD gives the rate of change of its phase along the flame front. Indeed, 
Eq. (PPD can be rewritten as 



dr, {v+f + 



(33) 



Then Eqs. (p2|), (^31) , (H), and ( p6|) allow one to express the right hand side of Eq. (|33|) 



in terms of the function /(r^), while Eq. (^) plays the role of the consistency condition 
which gives the equation for the function /(?]) itself. 
Specifically, Eqs. (0), (0), and (|2|) imply that 

= {vt? + Hf = ^{(^')' + i^')"} + (34) 

with some constant C. As in Eq. (|T9|) above, this constant is fixed up to the second 
order in ^ — 1, since to this order the flow is potential downstream, = Vp, uj~^ = uo^. 
The actual value of C can be found calculating Eq. (^) at some particular point at the 
flame front. At the tube walls, e.g., one has v~ = 1, = 6, ui^ = 0;+ = 0, therefore, 

C = e{e-1). (35) 

We will see in Sec. ^ that from the point of view of the large flame velocity expansion, 
another choice of the constant C is more appropriate: C = 6{6 — 1)V'^. This differs 
from that of Eq. ( |5BD only in the third order in 6' — 1, since V = 1 + 0{{6 — 1)^) (see 
Sec. |3]^ below). 

Next, solving Eqs. (|2^) , (|23|) with respect to v~^, u!~^, and using Eqs. (^4|), (^5|) , we 
derive the following expression for the phase derivative 

- _ pe^\ (r,- J. 

(36) 













e 


\V + v-y + {Hv-y + {6 - 1) 





Eqs. (p5D and (Pq) imply that the function v = v {rj) entering this equation can be 
represented in terms of f{ri) by formally inverting the operator (1 — f'H) : 

d- = {l-rH)-\N-V). (37) 

The right hand side of this equation can be represented as a formal series in powers of 
the operator f'H : 

(1 - fH)-\N -V) = {1 + f'H + f'H f'H + ---){N- V). 
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Now, to obtain the equation for /(r/), we have to use the remaining equation (29). 
In terms of the quantities Q and 0, it takes the form 

(Im + Re^) {ne'^ - OV) = 0, (38) 

where KeF and ImF denote the real and imaginary parts of the complex function F, 
respectively, and we used the fact that the Hilbert operator is real (see the Appendix). 
To combine Eqs. (|3^) , (|36D , and ( ^8l) into one, one has to solve the latter with respect 
to 0. This can be done as follows. 

Noting the relations ImiF = KeF, KeiF = —ImF, we can rewrite Eq. (|38D identi- 
cally as 

Im( 1 + iH) (fie*"^ -eV) = 0, (39) 



or as 



Re{i - H){ne''^ - 9V) = 0. (40) 

Acting by the operator H on Eq. ( ^0] ) from the left, and using the identity ( p^ of the 
Appendix, this equation can be rewritten also as 

Re{iH + 1) (fie^*^ ~ OV) = 0. (41) 

Together, Eqs. ( |39D and ( |4l| ) imply that 

{1 + iH){ne'^ - eV) = 0. (42) 



Let us consider the structure of Eq. (H2|) more closely. 



Lemma: All solutions of the equation 

(1 + zH)X = (43) 

span the ring. 

Proof: First, let us show that if X is a solution of Eq. (^), then its square X"^ also 
is. Indeed, using the identity (p8| ) of the Appendix, one has 



2iHX'^ = 2H{XHX) = {HXf - = -X^ - X^, 

therefore, 

{l + iH)X'^ = 0. 

In view of the linearity of the Hilbert operator, the sum Xi + X2 of two solutions Xi, 



X2 of Eq. (43) also is a solution, 

(l + zi^)(Xi+X2) =0, (44) 
and so are the squares Xf, X|, and (Xi +X2)^, as we just proved. We have, therefore, 
= (1 + iH){X, + X2f = (1 + iH){Xf + 2X1X2 + X|) = 2(1 + z^)(XiX2), 

or 

(l + 2i^)(XiX2) =0. (45) 
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In particular, it follows from Eq. (|45|) by induction that for any solution X of Eq. (|43|) 
and any positive integer n, 

{l+iH)X'' = 0, n&N. 



In mathematical terminology, the properties ( p4|) and (^51) reveal the ring structure of 
solutions of Eq. (|i3| ) . 

This result can be used to solve Eq. (|42|) with respect to 0. Namely, assuming that 
|f2e*'^/6'V^ — 1| < 1, and taking the infinite sum of the powers 



'W 



- 1 



multiplied by the factors -, we obtain 



1 ffle' 



- 1 



'l + iH)ln{l 



[l + ^H)[\n—+^^ 



0, 



or 



'W 



- 1 



This solution is then analytically continued for all values of fl, 



(46) 



(47) 



With the help of Eqs. ( p4D - (|36|), and (|4^) we can now write the equation for the 
function /(r/) we are looking for 



e d 



HI 



V 



9-1 
N 



{V + vf + {HvY + {9 -I] 



(4J 



with the denotation 



v^{l-fH)-\N-V). 



(49) 



Before we proceed to the investigation of this equation, let us consider the question 
of its accuracy. 



3.2 Accuracy assessment 



We derived Eq. ( ^81) using only exact transformations of the equations 
Therefore, its accuracy is determined entirely by the accuracy of the underlying model 
assumption expressed by Eq. (0), which can be estimated as follows. 

As we have mentioned above, the exact equation for the stationary flame front can. 



in principle, be obtained from the system of equations (^ - (p!l[) , of which only the 
?7-component of the Euler equations, Eq. (0), is not present in the set (*) of equations 
describing our model. The question of the model accuracy, therefore, is the question of 
the accuracy to which Eq. (|^ is satisfyed by the model solution. To answer the latter. 
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we need an explicit expression for the vorticity produced in the flame, which will be 
derived presently. 

It was shown in Ref. |jl2| that with the help of the conservation laws (H) - (|10D and 
the evolution equation (|llj), the value of the vorticity just behind the flame front can be 
explicitly expressed in terms of the fuel velocity. Namely, with the help of Eqs. (5.32) 
and (6.15) of the work [12|, the jump {cr}l of the vorticity 

du dw 
dr] 

across the flame front can be written, in the 2D stationary 

0-1 1 



= {dw. + f'Du. + -bf) , (50) 



where 



Differentiating the evolution equation (|TI]) and writing Eq. ( [5(]| ) longhand, expression 
in the brackets can be represented as a total derivative 

£)u>_ + /'£)«_ + Id/' s w'_w_ + + + /'«'_»_ + ^ 

= ^ in- - N)' ^ (N^ - iy_ ^ , (ul + (52) 

2 N N ' N 2 ^ ' 

Since the flow is potential upstream, we obtain the following expression for the vorticity 
just behind the flame front 

^+ = -^(«' +^-)'- (53) 

With the help of this equation, we can now show that Eq. (|^) is actually satisfied by 
the model solution with the accuracy of the sixth order in a = 6' — 1, as a ^ 0. 
First of all, the following estimates can be readily obtained [Q: 

f' = 0{a), u,wr.{f'f = 0{a'), 

and, more generally, 

'^ = 0(a)0(F), 

for any functional F = F[f{ri)], since the amplitude A of a perturbation of the flame 
front with the wavelength A is of the order aA |^, ^ ^. In particular, it follows from 
Eq. (H) that0 

a = O(a^). (54) 



^In the general 3D case, the same resuh follows from the formula 

(61-1)?;* r fi 
{rotv}± = - [n,Vtv*.], 

where [ , ] denotes the vector product, n the unit vector normal to the flame front (pointing to the 
burnt matter), v*_ the tangential to the flame front component of the velocity, and Vf differentiation 
in the direction v*l . 
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Next, we rewrite Eq. (|^) identically 

du dw 1 9n 

U— + W- ua = — — . (55) 

or] Of] p or] 

Setting ^ = /(?7) + in this equation and using Eq. (^3|), we obtain 

Using the above estimates and taking into account that Uj^ = 6V + u+, V = 1 + 0(a^), 
Eq. (|56D can be rewritten as 

+ +(g-l)# "-+"- =-^ (^] +0{a% (57) 



\dr] 2 J ^ dr] 2 

On the other hand, similar transformations of Eq. (^ give, with the same accuracy, 

/a«= + »n _ /anx ^^^^.^^ ^^^^ 



Finally, taking the sum of Eq. (^) and Eq. (|58|) multiplied by /', and noting that 

M^ + W;2 =y2_ +^2 ^^2 _ n+=n_-(0-l), 

we get 



^ l2 



ut + wi]+Il^ \ = 0{a'). (59) 



Since the flow is potential upstream, the left hand side of Eq. (|59D is zero. We conclude 
that Eq. (^), and therefore Eq. (0), is satisfied with the accuracy of the sixth order in 
a. 

It is worth to emphasize that the model equation (^) has not been used in the 
derivation of Eq. (|5PD. Therefore, the latter holds true whatever model is considered, 
provided that this model respects all the conservation laws at the flame front, as well as 
the flow equations up- and downstream. Furthermore, Eq. ([THp is mo del- independent 
up to the second order in a, since to this order the flow is potential downstream. 
Extended to all values of 6, Eq. (O) thus provides the simplest model satisfying the 



above-mentioned requirements. 

Finally, considered as an equation for the quantity {v?_ + w^)', Eq. (|59|) determines 
it with the accuracy of a^. On the other hand, since the left hand side of Eq. ( ^8]) is 
proportional to the same quantity, is the accuracy estimate for Eq. (|^) as well. 

It remains only to make the following three important remarks. 

Rl) The above accuracy estimate is obtained from the analysis of differential equa- 
tions governing the fluid dynamics. These equations should be complemented by ap- 
propriate boundary conditions. As we have already mentioned, the problem of flame 
propagation is essentially nonlocal; this non-locality shows itself in the fact that the 
boundary conditions for the burnt matter, together with the boundary conditions for 
the flame front itself, are invoked in the course of integration of Eq. (|48|) . By itself, this 
equation is independent of the boundary conditions, since it is obtained by the direct 
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substitution of Eq. (^) into the Euler equation (Q) for the burnt matter, written in 
the form (|36|). Thus, the fact that the consistency condition ( pQ]) can be resolved with 
respect to the phase is crucial for the above accuracy estimate. How the boundary 
conditions are actually taken into account in the course of integration of Eq. (|i8| ) will 
be shown in Sec. |4.1| . Closely connected to this is the remark 



R2) It was mentioned in Sec. |3]3 that the value (|35|) of the constant C entering 
Eq. ( p4D is fixed up to the second order, since to this order the flow is potential down- 
stream. One might think that the ambiguity in C at higher orders spoils the above 
O(a^) accuracy estimate for Eq. (^). It is not difficult to see, however, that within its 
accuracy, Eq. (|48|) is not affected by this ambiguity. Indeed, since the numerator in its 
right hand side is of the order a^, a third order change in C in the denominator gives 
rise only to terms of the order a®, and the same is true for the left hand side. This is 
what we should have expected, since as we proved above, Eq. (^S|) is the true equation 
with the accuracy a®, and as it is it must be independent of the model particularities 
in the higher-order completion of C. 

R3) As we mentioned in the points Rl), R2), equation (^8[) depends neither on the 
boundary conditions, nor on the higher order completion of the constant C. However, 



given the boundary conditions, C is ffxed upon integration of Eq. (|8). The point is 
that the boundary conditions for the ffame front together with the boundary condition 
for the ffuid velocity downstream imply more strong restriction on the value of C than 
that given by Eq. (|35|) to the second order in a. This is because the ffow structure near 
the ending points of the ffame front at the tube walls can be completely determined 
with the help of the boundary conditions. Thus, the initial choice of the constant C 
is effectively corrected in the course of integration of Eq. by the choice of the 
integration constants, appropriate to the given boundary conditions. In view of this 
fact, in practice, it is more convenient to work with Eq. (^Tf), the undifferentiated 
version of Eq. (^), from the very beginning, and choose the constant C from the 
requirement that the amplitude and phase take the boundary values prescribed by 
the boundary conditions. We follow this way in Sec. ^ below. Yet another example of 
determining the constant C is given in Sec. |[ 

3.3 Account of the finite flame thickness 



We will show in this section how the considerations of Sec. |3.1| can be generalized to 
take into account the effects due to finite fiame thickness. 

In the case of small but nonzero fiame thickness, the conservation laws at the fiame 



front read 112 



M+ - M_ - f\w+ - W-) = {6- l)N, (60) 
W+-W- + f'{u+ -uJ)=e\nd (dw^ + f'Du^ + ^^/') , (61) 





w'j" ^2DN -L_\ (62) 
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while the evolution equation 



N 



'9-ldri 



(63) 



where e is the small dimensionless ratio of the flame thickness to the characteristic 
length of the problem, and operator D is deflned in Eq. (pl|). 

As in Sec. one should first exclude f{r)) from the arguments of the flow variables 
entering these equations. Using the continuity equation (|^) and taking into account 
potentiality of the flow upstream, we write 



di] 
i)w. 




\drij_ \dri^ 



du 
dr] 



Df=(Du) - — 




■Df, 
Df 



(64) 
(65) 
(66) 



Introducing notation for the flow variables as in Sec. |3.1| , the above equations can be 
written in the form of the set (*), namely, up to the terms of the fourth order (which 
turns out be sufficient for expanding Eq. ( ^8| ) up to the fifth order). 



/ u+-ti- = {6-1)1 N \ 

W+-W- = -f'{e -1)/N + e\nef'f" 

n+-n- = -{e-i) + e{e-i)f" 

\u--f'w~ = N -eeine/{e -l){dw-/dri + du~/drif' + f") J 



(**) 



--m, 



all other equations remaining the same as in the set (*). Following the reasoning of 
Sec. all the fiow variables can now be considered ^-independent, which fact is 
expressed by the special designation (t;,ci;,7r) for the variables {u,w,Il), respectively. 

Thus, we see that in the case of nonzero fiame thickness, Eq. (|^) for the amplitude 
of the potential mode modifies to 



0{{v-)' + {uj~)'}-2ee{0-l)f" 



1). 



(67) 



The corresponding expression for the rate of change of the phase 0, which we do not 
write explicitly because of its complexity, can be obtained by substituting the equations 
(**) into Eq. (^). Differentiating Eq. (^Tj), and using the expression for 0' together 
with Eq. (^) for the amplitude, one obtains an equation for the function /(?]), which 
generalizes Eq. (^) for the case of nonzero flame thickness. 

The accuracy analysis for this equation is very complicated and will not be carried 
out here. It should be noted, however, that such an analysis is superfluous to a consid- 
erable extent. Indeed, the finite fiame thickness is mainly taken into account in order 
to provide a short wavelength cutoff for the spectrum of the flame front perturbations, 
which insures the existence of a stationary flame conflguration. On the other hand, to 
the leading order in a, the form of the e-corrections to the equation for f{ri) is known 
in advance from the linear theory of Pelce and Clavin 0. It will be shown in the 
next section that, along with the nonlinear ^-corrections, the equations obtained above 
correctly reproduce their result for the cutoff wavelength. 
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4 The small {9 — 1) expansion 



In its general form, Eq. ( ^8|) is very complicated. It is a highly nonlinear integro- 
differential equation, which can be solved only numerically. On the other hand, as we 
showed in Sec. |3.2| , this equation correctly approximates the exact equation for the 
flame front of zero thickness up to the fifth order in a. We now turn to expanding 
Eq. (|48|) , generalized to the case of nonzero flame thickness, up to this order. As the 
result, a much simpler equation will be obtained, which generalizes the Sivashinsky 
equation taking into account vorticity production in the flame. 



4.1 Fourth order equation for the flame front 

We will see below that to the order being considered, Eq. (^) can be integrated 
with respect to the phase, which implies that both sides of Eq. (^81) are actually full 
derivatives. However, it would be hasty to simply omit these derivatives: integration 
of Eq. (^81) requires careful account of the boundary conditions. Since this equation 
is integral, it is not sufficient to impose only one condition to fix the constants of 
integration in its left and right hand sides. Boundary conditions for the burnt matter 
and for the fiame front itself actually supply two independent conditions to be used to 
fix two arbitrary constants: one additive constant in the phase of the complex function 
(p^), and one multiplicative constant in its amplitude. 

We proved in Sec. |3.2| that within the accuracy of the sixth order in a, Eq. p8|) is 
independent of a particular completion of the constant (|35D beyond the second order. 
However, different choices of C correspond to different values of the integration con- 
stants in the phase and amplitude. It turns out that in the case of ideal tube walls, 
the particular choice ( pSj) made in Eq. ( ^8[ ) implies that integration of this equation 
gives exactly Eq. ( ^T)) where 0o = 0, 0o being the value of the phase at the tube 
walls. We will prove this statement only for zero-thickness flames. As we pointed out in 



Sec. |3.3| , nonzero flame thickness should be taken into account mainly in order to pro- 
vide the short wavelength cutoff for the flame front perturbations. On the other hand, 
e-dependence of the constants of integration is of little interest, and will be neglected 
in what follows. Let us now turn to the proof of the above statement. First of all, since 
the downstream flow is potential up to the second order, v"*" = v+ -|- O(a^), the phase 
(po = arg(Up -|- iWp) = 0{a^), since = at the walls. Furthermore, noting that 
u'^ + w'^ = (v_ ■ t)^ -|- (v_ • n)^, where t and n are the tangential and the normal to the 
flame front unit vectors, respectively, taking into account that v_ • tA^ = W- + f'u-, 
v_ ■ n = 1, and using Eq. (^) we obtain 

e-1 d {w. + fu.y 



2eN dri m 



It follows from this formula that a = at the walls, since the boundary conditions 
are assumed ideal. We conclude that the flow is potential near the walls, m+ = , 
= Wp, and thus 0o = 0. For the same reason, the amplitude f2 = 6^ at the walls. 



Now recall that the value (pSQ for the constant C was obtained in Sec. |3Tl| from the 
boundary conditions u~ = = 9, w~ = Wp = 0, valid up to the second order. 

We see, therefore, that in the case when the boundary condition for the flame front is 
/' = 0, Eq. (|35|) is extended to all orders in a, and thus Eq. (|47| ) is indeed the integral 
of Eq. p8|), satisfying the boundary conditions. 
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Let us now proceed with expanding Eq. (47) in powers of a. Since the numerator 
in the right hand side of Eq. (pSf ) is of the third order in a, it is sufficient to expand 
the denominator up to the second order: 



{(V^ + vf + [iivf - 2e{e - 1)/" + (^ - 1)} ' 



l-IW -1v 



11 



(6^ 



where W = V — 1 = 0{a^) is the flame velocity increase due to the front curvature, and 
the designation v~ is again reduced to v, for brevity. Substituting this into Eq. (|33D, 
we find 



h^ u/ (n f'^- 



l{i^„-/'(«-l)}^.0 



(69) 



It is not difficult to see that Eq. (|69|) can be integrated to give 



9 

Hv ,.,0-1 e\n9 „, „„ ^ ^ , 



9 



ON 



9 



(70) 



where 0o = for ideal tube walls, as we have shown above. One should remember, 
however, that 0o assumes different values under different boundary conditions. Thus, 
00 = 0{ea^), in general. It is worth to note also that this dependence of the flame 
equation on the boundary conditions for the burnt matter is a reflection of the essen- 
tial non-locality of the process of curved flame propagation, mentioned already in the 
Introduction and connected to the vorticity drift behind the flame front. 
Next, we expand the right hand side of Eq. ( ^T] ) 



n 

-H\n — 
9V 



-2e- 



-H{-\n 
I 2 



l + -(l^ + .) + ^-^ + 4^ 







■f" 



+ In 



1 + W 



(W + v) {Hvf - {W + vf 
~9 ^ 2 



—e- 



9 



(71) 



Substituting Eqs. ([7D|), ( |7TD into Eq. (^) and rearranging with the help of the 
identity (|98D gives 



9-1 



V 



W - H^) + '-{{9 - 1)/" + \n9Hifn} + O(a^). (72) 



On the other hand, one has from the last equation of the set {**)f\ 

(1 - fH)v = N-V- ^^{Hv' + v'f + /"), 



(73) 



^From now on, the terms O(a^) will be omitted, for simplicity. 
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which is the generahzation of Eq. (^) to the case of nonzero flame thickness. 
In the lowest, second, order in 6^ — 1, Eqs. ( |72D and (^) give, respectively, 

V = -^-^Hf, (74) 



y = \LL_W- ef". (75) 
The following equation for the flame front position, implied by Eqs. (0), (fTsP, 

^-^-W = -'^-^Hf + ef" (76) 

is nothing but the stationary part of the Sivashinsky equation [Q. 
In the next order, Eq. (^) becomes 

v=^-^{W-Hn + '-f". (77) 

Substituting this into Eq. ( [73D and expanding to the third order, we obtain the following 
equation 

Note that the expression ([79|) is just the first order approximation to the exact value 
of the cutoff wavelength 

given by the linear theory of the LD-instability [Q. 

Finally, at the fourth order, substituting Eq. (72) into Eq. ([73|) , and using the lower 
order equations (176]) and (|78D, we obtain 



= --^H!' + \ ("'""^ + 9 - l") /" + [m'f ) - S'H!"} . (81) 



We prefer not to expand the logarithm in the right hand side of this equation to 
make it transparent that the linear terms in Eq. (pi]) again correctly reproduce the 
corresponding terms of the linear theory. 
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4.2 Solution of the third order equation for the flame front 



The third order equation (|7q) is of the same functional structure as the Sivashinsky 
equation ([76|). Therefore, it can be solved analytically using the method of pole de- 



composition |I5|, |I6|. Considering the flame propagation in a tube of width b with ideal 
walls (for definiteness, the walls are taken to be the lines rj = and r] = b), we look for 
26-periodic solutions of the form 



fiv) 



2P 

a ^ In sin 

fc=i 



TT 



(821 



where the amplitude a and the complex poles 7]^, k = 1, ...,2P are to be determined 
substituting this anzats into Eq. (|78|). Since the function /(?]) is real for real rj, the 
poles come in conjugate pairs; P is the number of the pole pairs. Requiring the 2b- 
periodic solutions to be symmetric with respect to the reflection rj — >■ —rj, one can 
obtain periodic as well as non-periodic solutions to Eq. (^) in the domain r] G (0, b), 
satisfying the conditions /'(?7 = 0) = f'{ri = b) = 0. 
Using the formulae^ 



Hf = -7^ XI U + ^ sign(Im r]^) cot 



2P 



2b 



k=l 



TT 

2b 



sign(xj = 



X 



\x\ 



cot X cot y = —1 + cot(x — ?/)(cot y — cot x), 



^3) 



it is not difficult to verify that Eq. (|78D is satisfied by /(?]) taken in the form of Eq. (|82D 
provided that 




^(^ + 1) 2b \ 



2b 



^4) 



and the poles rjk, k = 1, 2P, satisfy the following set of equations 



2P 



i sign(Im rjk) + — J2 



m—1 



0, k = 1,...,2P. 



^5) 



It is seen from Eq. (^^ that for the tube width b > \^^\ the solution (^) is not unique: 
different solutions corresponding to different numbers P of poles are possible. To find 
the physical ones, the stability analysis is required which, of course, cannot be carried 
out in the framework of the stationary theory. However, as we have mentioned above, 
the functional structure of Eq. ( [78| ) is very similar to that of the stationary Sivashinsky 
equation (|76[) . Under assumption that the non-stationary versions of these equations 



are also similar, the stability analysis of Refs. [Ol^iHl will be carried over the present 
case. According to this analysis, for a given tube of sufficiently small width, there 
is only one (neutrally) stable solution. This solution corresponds to the number of 
poles that provides maximal fiame velocity. In addition to that, the pole structure 

^Since the application of pole decomposition to Eq. ( |78| ) is quite similar to that given in Refs. 
[T^ , we refer the reader to these works for more detail. 
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Figure 1: Dependence of the flame velocity increase on the inverse tube width, given 
by Eq. (B^) for the case 9 = 3, and the results of numerical experiments [5]. 



of the stable solution is such that the poles form a vertical alignment in the complex 
r^-plane, sharing the the same common real part. For such a "coalescent" solution, a 
simple upper bound on the number of poles can be obtained from Eq. ([85|). Namely, 
for k = ko with rj^^ uppermost, one has 



2P 



— — coth 
2b 



TT 



— (Im r]ko - Im r]„ 



> ^(2P- 1). 



m^kQ 



(The equality holds, if Im rjko = oo) Then it follows from Eq. ( ^If ) that the maximum 
of the flame velocity corresponds to the maximal number of the pole pairs 



Prr 



Int 



Ad) 



Int(x) denoting the integer part of x. Thus, the flame velocity increase Wg of the stable 
solution can be written as 



P A(^) 
' 2b 



P AW' 
2b 



(86) 



where 



46(6 + 1)' 



^7) 



Fig. compares the theoretically predicted dependence of the flame velocity increase 
on the inverse tube width, given by Eq. (^) for the case ^ = 3, with the results of 
numerical experiments [^. Dependence of the maximal flame velocity increase Wmax 
on the gas expansion coefficient, given by Eq. (^), is represented in Fig. |^. For 
comparison, we show also the corresponding dependence calculated with the help of 
the Sivashinsky equation. We see that the third order equation (^) gives reasonable 
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Figure 2: Dependence of the maximal flame velocity increase Wmax on the gas expansion 
coeflicient, given by Eq. (|87D (solid line), and calculated on the basis of the Sivashinsky 
equation (dashed line); the marks are according to the results of [5]. 

description of flames with ^ ~ 3, while for larger values of the expansion coefficient 
it overestimates the influence of vorticity, produced in the flame, on the flame front 
curvature. In the latter case, therefore, the more accurate fourth order equation (|ST| ) 
should be used instead of Eq. (^). Detailed investigation of Eq. (0) will be given 
elsewhere. 

5 Large flame velocity expansion 

As we have seen, Eq. (^) considerably simplifles in the case of small a. For sufficiently 
narrow tubes, it gives results which turn out to be in a reasonable agreement with the 
experiment already at the third order in a. Let us now consider the opposite case of 
very wide tubes, i.e., tubes of width large compared to the cutoff wavelength. As it 
will be shown presently, under a certain burning regime, Eq. (^) can be written in a 
much simpler form in this case too. 

Widely employed in modern jet engines is the process of the so-called fast flow 
burning (see, e.g., ||^, Ch.6, §1). This regime is characterized by a large stretch of the 
flame front along the tube, since the flame velocity is proportional to the flame front 
surface (in 2D case, to the front length). Indeed, equating the fuel flow at ^ = — oo 
to that through the flame front, using the evolution equation (0), and taking into 
account the boundary conditions w = 0, /' = 0, one has 
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Figure 3: Schematic representation of the flame front shape in the fast flow burning 
regime. 



or 

v = i»i. 

where / is the length of the flame front, n the unit vector normal to the flame front 
(pointing to the burnt matter), and h is the tube width. 

Under these circumstances, it is natural to develop an expansion of Eq. (|48|) in 
powers of the inverse flame velocity, l/V. This can be done as follows. 

Let us note, first of all, that large value of V implies, in general, that the quantity 
V (and, therefore, the quantity Hv - the 77-component of the fuel velocity at the flame 
front) is also large. Indeed, as it follows from the evolution equation (0), v = 1 — V ai 
the tube wall, in view of the boundary condition /' = 0. Now, we make the assumption 
that in the bulk, i.e., for all t] except for a small region near the points = 0, 6, the 
quantity v is actually of the order 0(1), rather than 0{y). Then the integral quantity 
Hv = 0(1) + 6/b 0{V) = o{V), since the size 6 of the region where t; ~ is assumed 
to be small compared with the tube width h, 5 <^h (see Fig. 

Expansion of Eq. ( ^8] ) in powers of 1/V is straightforward. Below, this expansion 
will be carried out only to the lowest nontrivial order. As in Sec.^, Eq. ( pO]) for the 
phase integrates at this order. Therefore, it is more convenient to work with the integral 
version (|47|) of the equation (^) from the very beginning. 

Let us first consider the question concerning the form of the quantities f2, entering 



this equation. As we saw in Sec. |3.1| , the value (^) of the constant C in the expression 
(P^) for the amplitude VL is fixed up to the second order in a. Furthermore, it was shown 
in Sec. |3.2| , that within the accuracy of our model, the form of Eq. ( ^8]) is independent 
of a particular completion of C beyond the order a^, and that the value of O, together 
with values of the arbitrary constants of integration of Eq. (|48|), can be deduced from 



the boundary conditions. In contrast, the boundary conditions cannot be used directly 
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in the framework of the 1/V^-expansion, since this expansion breaks down near the tube 
walls. Instead, the value of the constant C as well as the constant (po of integration 
of Eq. (|33|) will be determined from the requirement of self-consistency of the limiting 
transition ^ oo in Eq. (pTF]). 

Under assumption v, Hv ~ o{V), it follows from Eq. (^) that 0' ^ in the 
bulk when V ^ oo; therefore, = 0o = const. On the other hand, the amplitude 
Q y/OV^ + C in the same limit. Thus, Eq. (^Tj) gives in the lowest order in 1/V : 



^0 



-if In 



9V 

This equality is only consistent with the properties of the Hilbert operator if 
and therefore, 

00 = 0, c = e{e-i)v^. 



Note that in the case of small a, the values (pq) for the constants C, 0o coincide 
with those found in Sees. lO, 0|, respectively, up to the order a^. On the other 



hand, expansion of Eq. ( ^Tj) in powers of v/V has the strict validity in this case, since 
V ~ O(a^), V = 1 + O(a^) for all t] G (0, b). We conclude, therefore, that the first 
order of the large V expansion of Eq. ( ^TD must automatically reproduce the second 
order of the small a expansion of this equation, i.e., the Sivashinsky equation ([76|). 
Substituting Eq. (|^) into Eq. (|47| ) and expanding to the first order in v/V gives 

V {Hv - /'^ 



-H— 

e^v^ ev 

or 




V = — -H 4? • (89) 



Since the flame is highly stretched in the fast flow burning regime, the slope /' of the 
flame front is large. However, it follows from Eq. (^) that the quantities v and Hv 
remain of the order 0(1), since f'/N < 1. Our initial assumption is thus confirmed. 
Substituting Eq. (^) into the evolution equation (|^), we obtain 

We see that in the case of small 9 — 1, Eq. ( ]90| ) does reproduce the stationary part of 
the Sivashinsky equation (^) in the case e = 0. 

^The operator (^) is badly defined for f{ri) = C = const, except for C = 0. However, one can 
formally extend this definition to include nonzero values of C by setting HC = (3C with some number 
(3 which must be one and the same for all C, in order to preserve the linearity of the Hilbert operator. 
Furthermore, the properties (|9^) and (|9^) are preserved only if /3 = ±i. Since the quantities 0, are 
real by definition, the choice (|88|) is unique. 
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Finally, it is not difficult to take into account the influence of the small flame 
thickness. First of all, we note from Eq. (p9| ) that although the quantity v = 0(1), its 
r]-derivative 




since f\N ^ f ~ V. It follows then from Eqs. (H-dl) that Du^, Dw- = 0{1/V). 
Resolving Eqs. (0), (|6T|) with respect to w+, we see that the 5-corrections to 
these quantities are only of the relative order 1/V'^. Analogously, the e-correction in 
the pressure jump (|6^ ) is 0(1), which implies the same 0(1/V^^) relative correction in 
the amplitude Q. Thus, we conclude that Eq. (|47|) remains unchanged in the ffist order 
in 1/V, and so does, therefore, its consequence, Eq. (pOf). Substituting the latter into 
evolution equation (|63D , we obtain 

It is interesting to note that the term describing the influence of finite flame thickness 
turns out to be linear in /. 



6 Discussion and conclusions 

We have shown that in the stationary case, the asymptotic expansion of the nonlinear 
equation for the flame front position can be pushed beyond the second order in ^ — 1, 
at which the gas flow is potential on both sides of the flame front, to take into account 
vorticity drift behind the flame front. This expansion has been carried out explicitly; 
for the case of ideal tube walls, it is given by Eqs. ([T8|), (^Tj) at the third and fourth 



orders, respectively. Remarkably, the third order equation, which describes influence 
of the vorticity on the flame front structure to the lowest nontrivial order, turns out 
to have the same functional structure as the Sivashinsky equation. As we showed in 



Sec. |4.2| , it gives results in a reasonable agreement with the numerical experiments on 
the flame propagation in tubes for the case of flames with ^ < 3. 

It should be mentioned that the nonlinear equation for the flame front, derived 



in Ref. [0, also gives satisfactory description of the stationary flame propagation in 
narrow tubes. We showed already in the Introduction that the approach of Ref. [|l3 
is self-contradictory. Still, one might imagine that despite an erroneous derivation, the 
resulting equation is correct. However, the results of Sec. f4.1| show that in the case 
6^1, asymptotic expansion of the true nonlinear equation is quite different from that 
of the equation proposed in Ref. [|1^]. The latter is incorrect, therefore, already in the 
case of small gas expansion. 

From a more general point of view, the theory of flame propagation in the fully 
developed nonlinear regime cannot be formulated in the way the Sivashinsky equation 
1^ or the Frankel equation |^ are formulated. This is because the assumption of 
potentiality of the flow downstream, employed in these works, renders the relations 
between the flow variables local, allowing thereby the formulation of equation for the 
flame front position in terms of this position alone. In contrast, account of the vorticity 
production in the flame makes the problem essentially nonlocal, since, e.g., the value 
of the pressure field on the flame front is a functional of the velocity field in the bulk, 
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which imphes, in particular, that the equation for the flame front position must be a 
functional of the boundary conditions. Indeed, as we saw in Sec. |4.1j , the boundary 
conditions on the tube walls are invoked in the course of derivation of the asymptotic 
expansion of this equation already at the third order. 



Finally, with the help of the general equation (0), highly nonlinear regimes of the 
stationary flame propagation can be considered. We would like to remind that this 
equation respects all the conservation laws to be satisfied across the flame front, by 
construction. Thus, despite the fact that the vorticity drift behind the flame front 
is taken into account in this equation on the basis of the model assumption (pTS]), 
unjustified for arbitrary 6*, one may hope that it gives at least qualitative description. 
We showed in Sec. |] that in the particular case of the fast flow burning, equation (|i8| ) 
(as well as its generalization to the case of nonzero flame thickness) can be highly 
simplified by expanding it in powers of the inverse flame velocity. The result of this 



expansion to the first order is given by Eq. (|9T|). 
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7 Appendix 

For the sake of completeness, we give here a brief account of the properties of the 
Hilbert operator, used in the text. 

Given a sufficiently smooth integrable function /(t]), t] G (— oo,+oo), the Hilbert 
operator H is defined by 

(HfM = -p.v. / dcP^, (92) 

— oo 

"p.v." denoting the principal value. By definition, operator H is linear, i.e., for any 
complex numbers Ci, C2, 

and real, i.e., 

{HfY = Hf\ 

where F* denotes the complex conjugate of F. 

It is convenient to introduce the usual scalar product of two functions fiirf) and 

+00 

(/i,/2)= J dr^mv)f2{v)- 

—00 

Then, changing the order of integration, one has 

{f„Hf,) = - / rfr^/r(r])p.v. / dCp^ = — / rfC/2(C)p.v. / dr^l^ = -{H hJ,), 

TV J J C — ri 71 J J 71 — C 
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or 



(/l,^/2) = -Wl,/2), (93) 

i.e., the Hilbert operator is anti-Hermitian. 
To prove the operator identity 

= -1, (94) 

it is convenient to represent the right hand side of Eq. (^) as the integral over the 
contour Ci = U in the complex ?7-plane (see Fig. H) 

Then the square of the Hilbert operator takes the form 

(H^f)i,)^l.lJ^li.m, (96) 

47r^ J z — 1] J z — z 

where the contour C2 = C2 U C2 comprises Ci. Changing the order of integration in 
Eq. (|96D , using the formula 



/ a 



dz 1 ^ z — ti 

m ■ 



{z — ri){z — z) z — 7] z — z^ 

and taking into account that the logarithm gives rise to a nonzero contribution only if 
the arguments of the functions z — r] and z — z run in opposite directions when z runs 
the contours Cf , we obtain 



1 t a r^j(^) _ 1 



4:71^ J Z — ri J Z — Z 4:71^ 

Ci ' C2 

= ^fdz-^^'^ 



I dzl^ - 2.^ [ dz^^ 

J Z — 7] J Z — 7] 

V c+ ) 



c 



It^ 1^ z — rj 

and therefore, 



{H'fM = -fiv). 

Next, let us consider the quantity 

1 r.^ns) 



z — ri J z — z 

Ci ' C2 



,j,mu,m^2^f,ui^f,,m. (97) 



47r^ J z — T] J z — rj 47r^ J z — rj J z — z 

The first term in the right hand side of Eq. ( pTl) is just {{Hf ){r])}'^. The second term 
can be transformed as follows. The contour Ci of the 5-integral can be moved down 
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Figure 4: Contours of integration in Eqs. (|95D-(p7D. 

pass the contour C^, the pole at z = z giving rise to the extra term —2TTif{z), z G C^. 
Likewise, Ci can be moved up pass C2, with the extra term 2nif{z), z G Thus, 

471"^ J Z — 7] J Z — Z A-K^ J Z — 7] J Z — Z 

C2 Ci Ci C2 



2 



\P2" c+ J 



1 /,./(£l 

At:^ J z — ri J z — z 

Ci ' C2 



f dzl^ = {H{fHf}){r^) + f{r^). 
Jc z — n 



2tt Jc z — 1] 

Substituting this into Eq. (0), we arrive at the following identity 

2H{fHf} = {Hff - f. (98) 

Finally, let us establish the connection 

<! = -dH (99) 

between the 2D Landau-Darrieus operator $ and the Hilbert operator. The former is 
defined by 

+00 

(l>/)(r/) = / dk\k\Ue^''\ 



"sign" 



fk being the Fourier transform of /. Using the Fourier representation of the function 

1 \ 



1 / 1 

—00 ^ 

one has 

+00 

($/)(r/)= J dkksign{k)fke"'^ 



-iO C + iO 
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— oo 
hoo 



T^p.v. / dC- 



c 



— oo ^ — oo 

or 

Eq. (^) is convenient in estimating the orders of magnitude of expressions involving 
the operator $. In view of Eqs. (01) and (CT), one has 



{Hf,Hf) = ~{H'f,f) = {f,f), 
i.e., the Hilbert operator has unit norm, and therefore 

^f = 0{df). 
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